SELF-SIMILAR SOLUTIONS OF THE NON-STRICTLY HYPERBOLIC 
WHITHAM EQUATIONS FOR THE KDV HIERARCHY 

V. U. PIERCE AND FEI-RAN TIAN 

Abstract. We study the Whitham equations for all the higher order KdV equations. The Whitham 

o " 

£NJ ' equations are neither strictly hyperbolic nor genuinely nonlinear. We are interested in the solution of 

S-H 

O !■ the Whitham equations when the initial values are given by a step function. 

< 

1. Introduction 
| It is known that the solution of the KdV equation 



S3 



1.1) u t + 6uu x + e u xxx = 



> . 

has a weak limit as e — > while the initial values 

in 

u(x, 0;e) = u (x) 

1: „„ 



This weak limit is described by hyperbolic equations. It satisfies the Burgers equation 

(1.2) u t + (3u 2 ) x = 

until its solution develops shocks. Immediately after shock, the weak limit is governed by the 
Whitham equations EJ dH [15] 

(1.3) u it + Xi(u 1 ,u 2 ,u 3 )u ix = , i = 1, 2, 3, 

where the A«'s are given by formulae (12.121) . Equations ( 11.3j) form a 3 x 3 system of hyperbolic 

equations. After the breaking of the solution of (11. 3p . the weak limit is described by a 5 x 5 systems 

of hyperbolic equations similar to (11.31) . Similarly, after the solution of the 5x5 system breaks down, 

l 
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the weak limit is characterized by a 7 x 7 system of hyperbolic equations. In other words, for general 
initial data uo(x), one constructs the weak limit by patching together solutions of (11.21) . (11. 3ft . 5x5, 
7x7, etc systems in the x-t plane. 

The KdV equation (II. ip is just the first of an infinite sequence of equations. All these so-called 
higher order KdV equations can be cast in the Hamiltonian form 

du d 8H m 

where H m 's form a sequence of conserved functionals for the KdV equation. The small dispersive 
parameter e is hidden in H m . In particular, when m = 1, (11.41) is the KdV equation. 

The solution of the higher order KdV equation (11.41) also has a weak limit as e — > 0. As in the 
KdV case, this weak limit satisfies the Burgers type equation 

(1.5) u t + (-^j—u m+1 ) x = , 

m + 1 

where C m is given in (I2.20p . until the solution of ( 11. 5ft forms a shock. After the breaking of the 
solution of ( 11.5ft , the limit is governed by equations similar to ( 11. 3ft . namely, 

(1.6) u it + fi ( f l \ui,U2,u 3 )u ix = , i = 1,2,3, 

where /4"^'s are given in (12.151) . They will also be called the Whitham equations. As in the KdV case, 
after the solution of (II. 6p breaks down, the weak limit is described by a 5 x 5 system of hyperbolic 
equations. 

In this paper, we are interested in the solution of the Whitham equations for the higher order KdV 
( 11.41) with a step-like initial function 



1.7) u (x) 



For such an initial function, the solution of the Burgers type equation (II. 5p has already developed a 
shock at the initial time, t — 0. Hence, immediately after t = 0, the Whitham equations (11.61) kick 
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in. Solutions of (ll.6p occupy some domains of the space-time while solutions of (11. 5p occupy other 
domains. These solutions are matched on the boundaries of the domains. 

Equations ( 11.2ft and equations (11.51) are prototypes in the theory of hyperbolic conservation laws 
[TJ. Their solutions will generally develop shocks in finite times. The solutions can be extended 
beyond the singularities as the entropy solutions. 

Solutions of equations (11.21) or equations (11.51) . in the theory of the zero dispersion limit, are not 
extended as weak or entropy solutions after the formation of singularities. Instead, they are extended 
to match the Whitham solutions of ( 11.31) or ( 11. 6ft . For initial data (11.71) . the resulting solutions of 
the Whitham equations (11.61) will be seen to be more complex than those of (jl.3p in the KdV case. 

The KdV case with the step-like initial data (11.71) was first studied by Gurevich and Pitaevskii 
[3]. They found that it was enough to use the Burgers solution of (ll.2p and the Whitham solution of 
(11.31) to cover the whole x-t plane, without going to the 5 x 5 or 7 x 7 system. Namely, the space-time 



is divided into three parts 

(1) ~<-6, (2) "6<-<4, (3) - t >4. 
The solution of (11. 2p occupies the first and third parts, 

(1.8) u(x, t) = 1 when | < — 6 , u(x,t) = when f > 4 . 
The Whitham solution of (jl.3p lives in the second part, 

(1.9) Ul (x, t) = 1 , | = A 2 (l, us, 0) , u 3 (x, t) = 0, 
when — 6 < x/t < 4. 

Whether the second equation of (11.91) can be inverted to give w 2 as a function of the self-similarity 
variable x/t hinges on whether 

^(1,^2,0)^0. 

OUo 
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Indeed, Levermore [8] has proved the genuine nonlinearity of the Whitham equations (11.31) . i.e., 

d\ 

(1.10) -^( Ul ,u 2 ,u 3 ) > 0, 2 = 1,2,3, 

OUi 

for u\ > u 2 > u 3 . 

For the higher order KdV (jl.4p . equations (11.61) . in general, are not genuinely nonlinear, i.e., a 
property like (11.101) is not available. Hence, solutions like (jl.8p and (11.91) need to be modified. 

Our construction of solutions of the Whitham equation (ll.6p makes use of the non-strict hyperbol- 
icity of the equations. For KdV, it is known that the Whitham equations (11.31) are strictly hyperbolic, 
namely: 

Ai(ui,u 2 ,u 3 ) > A 2 (ui,u 2 ,u 3 ) > X 3 (ui,u 2 ,u 3 ) 

for u\ > u 2 > U3 [8J. For the higher order KdV (11.41) . different eigenspeeds of (11.61) . fi[ m \ui, u 2 , «3)'s, 
may coalesce in the region u\ > u 2 > u% [9]. 

For the higher order KdV with step-like initial function (11.71) . the space time is divided into four 
regions (see Figure 1.) 

(l)f<«, (2) a<|</?, (3)/5<^<4 m , (4)^>4 m , 

where a and j3 are some constants. In the first and fourth regions, the solution of (11.51) governs the 
evolution: 

u(x,t) = 1 where x/t < a and u(x, t) = where x/t > 4 m . 
The Whitham solution of ( 11. 6p lives in the second and third regions; namely: 

(1.11) m(x,t) = l, - = ^ { 2 m \l,u 2 ,u z ) , I = iif l \l,u 2 ,u z ) , 
when a < x/t < /3, and 

(1.12) Ul (x,t) = l, j =/i 2 m) (l,u 2 ,0) , u 3 (x,t)=0, 
when (3 < x/t < 4 m . 
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Equations (11. lip yield 

/4 '{l,u 2 ,u 3 ) = fM 3 '{l,u 2 ,u 3 ) 

on a curve in the region < u 3 < u 2 < 1- This implies the non-strict hyperbolicity of the Whitham 
equations (jl.6p for the KdV hierarchy. 

The m = 2 case has been studied in [9]. There, inequalities 

a M o M _ (m) q (m) 

(i.i3) ^5_<££ — ^ < cn^_ forUi>U2>M3>0 

tfw 3 2 « 2 — u 3 ou 2 

have played a crucial role in verifying that equations ( 11. Ill) or ( j 1 . 121) can indeed be solved to give 
the solution of the Whitham equations (11.61) when m = 2. 

For m > 2, inequalities (I1.13P are not valid any more. We therefore must use a different approach 
to solve the problem. The calculations are considerably more difficult than in the m = 2 case. This 
is mainly because q of ( I2.25P is a polynomial of degree m when m > 2 while it is only a quadratic 
polynomial when m — 2. 

The organization of the paper is as follows. In Section 2, we will study the eigenspeeds, /4"^'s, of 
the Whitham equations (11.61) . In Section 3, we will construct the self-similar solution of the Whitham 
equations for the initial function (11. 7p . In Section 4, we will use the self-similar solution of Section 3 
to construct the minimizer of a variational problem for the zero dispersion limit of the KdV hierarchy. 

In a subsequent publication, we will study the Whitham solutions for all the other step-like initial 
data. 

2. The Whitham Equations 

In this section we define the eigenspeeds of the Whitham equations for both the KdV ( 11.11) and 
higher order KdV (ll.4p . We first introduce the polynomials of £ for n — 0, 1, 2, . . . [Tj HJ [12] : 

(2.1) P n (C, ui, u 2 , u 3 ) = C +1 + a n ,iC + ■■■ + a„,„+i , 
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where the coefficients, a nj i, a nj2 , . . . , a njTl+ i are uniquely determined by the two conditions 

(2.2) - Pn^U^U, ) = ^ + 0(r3/2) ^ ^ 

and 

(2.3) r — ^ = o . 

■A* V (f _ n (f _ ^2) (£ - u 3 ) 
Here the sign of the square root is given by y(| — tti) (f — m 2 ) (| — ^3) > for £ > ui and the branch 
cuts are along (— 00, m 3 ) and (w 2 ,wi). 
In particular, 

(2-4) P (^,u 1 ,u 2 ,u 3 ) = £ + ao,i , Pi(^,Mi,u 2 ,m 3 ) = £ 2 - + u 2 + u 3 )£ + a h2 , 



where 



Here 



K{s) 

Ol,2 = 7i{ u l u 2 + «1«3 + M2«3j + 7(^1 + «2 + M3jao,l 
O 



«2 — «3 
S = 



Ui - U 3 

and K(s) and £?(a) are complete elliptic integrals of the first and second kind. 

K(s) and E(s) have some well-known properties [T0l[lT]. They have the expansions 



7T M s 9 2 ,l-3---(2n- 1). 2 _ 
- 1 + - + — s 2 + ■■■+ — _ . v _ '- V 



(2.5) K(s) = -[1 + - + — s 2 + ■■■ + (— ^ — + •••]. 

v; v; 2 L 4 64 v 2 • 4 • ■ • 2n ' J 

(2.6) £(a) = *[l- S --^s 2 l_ ( l-3---(2n-l) 

v ; v ; 2 L 4 64 2n - l v 2 • 4 • ■ • 2n ' 

for |a| < 1. They also have the asymptotics 

1 16 

(2.7) K(s) « -log 



2 ° 1- a ' 

1 16 

2.8) £(a) « i + -(i_ s )[lo gl _-l] 
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as s is close to 1. Furthermore, 

dK{s) E(s) - (1 - s)K{s) 



(2.9) 
(2.10) 



ds ' 2s(l-s) 
dE(s) E(s) - K(s) 



ds 2s 
It immediately follows from (12.51) and (12. 6p that 

The eigenspeeds of the Whitham equations (II. 3p are defined in terms of P and Pi of (12. 4p . 

X i {u 1 ,u 2 , u 3 ) = 12— r , i = 1, 2, 3 , 

Fo{Ui,u 1 ,u 2 ,u 3 ) 



which give 



K(s 

Xi(ui, u 2 , u 3 ) = 2(ui + u 2 + u 3 ) + 4(«i - u 2 ) 



(2.12) A 2 (mi, u 2 , u 3 ) = 2(«i + u 2 + w 3 ) + 4(w 2 - «i) 

A3(«i, «2, M3) = 2(«i + u 2 + M3) + 4(m 2 - M3) 



E(s) - (1 - s)lf(s) 



£(s) - K{s) 

In view of (I2.5ti2.8p . we find that Ai, A 2 and A3 have behavior: 
(1) At u 2 = u 3 : 



(2.13) 

(2) At ui = u 2 
(2.14) 



\x{ui,u 2 ,u 3 ) = 6«i , 

A 2 («i, u 2 , u 3 ) = X 3 (ui,u 2 , u 3 ) = I2u 3 - 6«i . 



Xi(ui, u 2 , u 3 ) = X 2 (m, u 2 , u 3 ) = Aui + 2w 3 , 
X 3 (ux,u 2 ,u 3 ) = 6u 3 . 
The eigenspeeds of the Whitham equations (II. 6p are 



(2.15) n\ '(u 1 ,u 2 ,u 3 )=A (2m +1)— -, 2 = 1,2,3 

Po{v>i,u>i,u 2 ,u 3 ) 
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The polynomial 4 m (2m + l)P m (£, w l5 u 2 , u 3 ) can be expressed as [2] 

(2.16) 4 m (2m + i)P m (£, mi, w 2 , w 3 ) = 2(f - - « 2 )(£ - «3)$(£, «i, Ua, «a) + Q(f , «i, «2, «s) • 
The function $(£, «) satisfies the boundary value problem for the Euler-Poisson-Darboux equations 



(2.17) 2(ui-Uj) 



duiduj dui duj 



. <9 2 $ 9$ 9$ 

< 2 - 18 > 2 «-"->^- = wr 2 ^' 

(2.19) *(«,«,«,«) = |^[C m « m ], 

where 

2 2m+1 2 m (2m+l)!! 
2.20 C m = = — i ;— ^- . 

Jo v 7 ! 3 * 

The function Q(£,u) is a quadratic polynomial in £; 

(2.21) Q(f,Ui,u 2 ,U8) = 2(^-U2)(?-W3) a + 2(f-Ui)(f-u 3 ) ~ 

+2(4 - UlM - U V n H ?(«l,tt2,«3)"o(c,ttl,«2,«3) 

and g(w) is the solution of the boundary value problem for another version of the Euler-Poisson- 
Darboux equations 

. d 2 q dq dq 

(2 - 22) ^-^adr. = £-w-> »''' = 1 ' 2 ' 3 ' 

L/ Li' ^ L/ LL j !_/ LL^ LL j 

(2.23) g(u,u,u) = C m u m . 

The solution of equations (l2.17H2.i9l) and that of (12.221) and (12.231) can be solved explicitly [2] . In 
particular, the solution of flZggp and f[2T23l is [ID] 



In 0/l^ / \ °w / / V 2 2 U1 ^ 2 2 " 2 ^ 2 j j 

(2.24) q(u 1 ,u 2 ,u 3 ) = — -= \ / === d/xdi/ . 



2^ 7_! 7_i ^(1-^(1-^) 
The speeds //^ 's of (I2.15P for m > 1 are connected to /4 s for m = 1, which are also given by 

(TO . 
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Lemma 2.1. For i = 1, 2, 3, 



'o on "" ! / 1 w v of , , m 9q(u 1 ,u 2 ,u 3 ) 

(2.25) /^ '(1*1,1*2,1*3) = -[A^U!, h 2 , w 3 ) - 2(i*i + w 2 + i* 3 )] 7^7 \- q{u 1 ,u 2 ,u 3 ) . 

Proof. We use (12151) . (12161) and (12T2T1) to write 

(m)/ \ Q(i*i, 1*1,1*2, 1*3) 

A*a («i, 1*2,1/3) = -57 7 

^0(1*1,1*1,1*2,1*3) 

, ooc x 2(i*i - i* 2 )(i*i - 1*3) dq(u 1 ,u 2 ,u 3 ) 

(2.26) = — — r + q{u x , 1*2, u 3 ) . 

2-0(1*1,1*1,1*2,1*3) wi*i 

In particular, when m = 1, since the corresponding q = 2(ui + u 2 + u 3 ), we obtain 

Ai(i*i, i* 2 , i* 3 ) = — ^7 s h 2(i*i + u 2 + u 3 ) . 

^(1*1,1*1,1*2,1*3) 

This together with (I2.26P proves formula (I2.25P for i — 1. The cases for z = 2, 3 can be shown in the 
same way. □ 

Lemma 2.2. [12] 

i. 

(m) 

(2.27) = - ^ , i, j = 1, 2, 3;i ^ j . 

<9 / P m (£, l*i, 1*2, 1*3) \ /tj m) (l*l,lt 2 ,l*3) ^ / -Fq(^, 1*1, 1*2,1*3) 



2.28 



^ VV(e-«i)(e-«2)(e-«3)y 4™(2m+l) ft* \y ( £_ Ul )(f _ U2 )(f _ U3 )_ 
/or z = 1, 2, 3 and r\ 7^ u 1( u 2 , u 3 . 

The following calculations are useful in the subsequent sections. 
Using formula (I2.25P for /* 2 and fi 3 and formulae (I2.12p for X 2 and A3, we obtain 

(2.29) /i 2 '(1*1,^2,1*3) -/*3 (1^1,^2,1*3) = ^ _ - (1 - g^] M ( Ml ' M2 ^ 3 ) > 
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where 

(2.30) M( Ul ,u 2 ,u 3 ) = [£ + (l - .)£]E - (1 - .)(£ + . 

We then use flZSD, (12101) and (12T22I) to calculate 

/ 231 n dM(u 1 ,u 2 ,u 3 ) 1 pi(u u u 2 ,u 3 ) 

du 2 2 u\ — u 3 

dM(u u u 2 ,u 3 ) 1 p 2 (u h u 2 ,u 3 ) 3 M(mi,m 2 ,m 3 ) 
i 2 - 62 ) a = o [h-{l-s)K + - , 

OU 3 2 U\ — U 3 I U\ — u 3 

where 

.9 d 

(2.33) pi(«i,it2,ti 3 ) = 2(ni-M 2 )^-div(g)-div(g) , P2(u 1 ,u 2 ,u 3 ) = 2(«i -w 3 )— div(g) -div(g) 

We next consider 

(2.34) F(tti,tt2,tt 3 ) := • 

Using formula (I2.25P for fi 2 and fi 3 and formulae (12.121) for X 2 and A3, we obtain 

F = 2 (1 ~ S)ir ^ I 2 * 9g 



E — (1 — s)if <9w 2 K-Edu 3 
, , _ v 2[ g (l-s)if ,_gg 

where we have used equations f!2.22j) in the last equality. Finally, we use the expansions (l2.5H2.6p for 
K and E to obtain 

(2.35) F(u h u 2 ,u 3 ) = -4[(2 - I a + ■ ■ - - «s) + (-7 + 0(s 2 ))p-] . 

4 OU 2 OU 3 4 (TO 3 

3. Self-similar Solutions 

In this section, we construct the self-similar solution of the Whitham equations (jl.6p when m > 2 
for the initial function (11.71) . The m = 2 result has already been obtained in {9j. Even in the m = 2 
case, the key calculations presented here are different from those in [9]. 
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Theorem 3.1. (see Figure 1.) For the step-like initial data Uq(x) of ( [i. 7\ ), the solution of the 
Whitham equations $1.6\) is given by 

(3.1) ui = l, x = /4 m) (l,u 2 ,u 3 ) t , x = /i^ (1, u 2 , u 3 ) t 
for at < x < fit and by 

(3.2) ui = l, x = fi ( 2 m \l,u 2 ,0) t , u 3 = 

for j3t < x < jt, where a = ^\l,u\u*), (3 = fx ( 2 m \l, u**, 0) and 7 = ^ m) (l, 1, 0) = g(l,l,0) = 4 m . 
i/ere ; is uniquely determined by the equation 

(3.3) p 1 (l,u*,u*) = 0, 
and zs uniquely given by the equation 

(3.4) /4 m) (l l «-,0)-A*S m) (l,«-,0) = 0. 

Outside the region at < x < 4 m t, the solution of the Burgers type equation U.5\) is given by 

(3.5) u = 1 x < at 
and 

(3.6) u = x > A m t . 

The boundaries x = at and 2 = 4 m t are called the trailing and leading edges, respectively. They 
separate the solutions of the Whitham equations and Burgers type equations. The Whitham solution 
matches the Burgers type solution in the following fashion (see Figure 1.): 

(3.7) Ui = the Burgers type solution defined outside the region , 

(3.8) u 2 = u 3 , 
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U1 








U2. 






X 








a B 


4 m t 



Figure 1. Self-Similar solution of the Whitham equations for m > 2. The curve 
defines the function x( u )- 

at the trailing edge; 



(3.9) 
(3.10) 



Ml = M 2 , 

u 3 = the Burgers type solution defined outside the region 



at the leading edge. 

The proof of Theorem 13.11 is based on a series of lemmas. 

We first show that the solution defined by either formulae (13. ip or (13. 2p indeed satisfies the 
Whitham equations ([HD [HO [13]. 

Lemma 3.2. (1) The functions u\, U2 and M3 determined by equations [3. 1\) give a solution of 
the Whitham equations $1.6}) as long as M2 and M3 can be solved from $3. II) as functions of x 
and t. 

(2) The functions U\, U2 and M3 determined by equations A3. 2\) give a solution of the Whitham 
equations U.6]) as long as M2 can be solved from Ii3.2\) as a function of x and t. 



Proof. (1) Mi obviously satisfies the first equation of (11.61) . To verify the second and third equations, 
we observe that 



(3.11) 



d^ 2 m) _ <9/4 m) 



du-) 
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on the solution of (13.11) . To see this, we use (I2.27P to calculate 



(m) 8X2 



^f^2 9« 3 / (m) (in) 



du 3 X 2 - A 3 
The second part of (13. lip can be shown in the same way. 

We then calculate the partial derivatives of the second equation of (13.11) with respect to x and t. 

1 = -5 tu 2x , = — tu 2t + /4 , 

which give the second equation of (11.61) . 

The third equation of (11.61) can be verified in the same way. 
(2) The second part of Lemma 3.2 can easily be proved. 

□ 

We now determine the trailing edge. Eliminating x and t from the last two equations of ( 13.11) 
yields 

(3.12) fx { r\l,u 2 ,u 3 ) - fxt\hu 2 ,u 3 ) = . 
Since it degenerates at u 2 = u 3 , we replace (13.121) by 

(3.13) F(l, U2, u 3 ) := ^\^u 3 )-^\l^u 3 ) = Q 

u 2 - u 3 

Here, the function F is also defined in (I2.34p . 

Therefore, at the trailing edge where U2 = u 3 , i.e., s = 0, equation ( I3.13p . in view of the expansion 
( I2.35p . becomes 

d 2 q(l,u 2 ,u 2 ) dq(l,u 2 ,u 2 ) 

8(1 - u 2 ) — ^ — 5 3 = . 

ou 2 ou 3 ou 3 

Since = and l^f = 3 ^^ on u 2 = u 3 because of (12.221) . this equation is exactly equation 

Lemma 3.3. Equation £) = has a simple zero, denoted by u* , in the region < £ < 1, 

counting multiplicities. Furthermore, £) positive when t; < u* and negative when £ > u*. 
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Proof. We first simplify the polynomial p\ of (I2.33p . In view of formula (12.241) for q, we use the fact 

that q is symmetric in m, u 2 and 113 to obtain 

,ai A \ ( dq(u u u 2 ,u 3 ) dq(u u u 2 ,u 3 ) dq(u u u 2 ,u 3 )\ 

(3.14) ^ — + — + — J \ U2=U3= , = C m U^ Ul ) , 



where 



1 m (i+M^ + l^£ Ml 



m—1 



(3.15) <7 °«--> = 57s7_ r w - 

We can then write p\ as 

(3.16) = c m [(i-o du f^ 1] - tf„(e, i)] • 

Denoting the function in the parenthesis of (I3.16P by p(£), we claim that 
(317) ^P(0) >0 ^Ml) <0 

for A; = 0,1, 2, ••• ,m-2. 
Obviously, 

Since, d ^f' 1 ' 1 is a positive function, this proves the second inequality of ( I3.17p . 
To prove the first inequality of (13.171) . we use formula ( 13. 15ft to calculate 

d k U (^l) , _ m(m-l)--(m-fc) f 1 „w-fc-§n , , A k Hll 
— — k=o- 2™+| ]_} M) 2 (l + /i) ^- 

The integral on the right can be evaluated using an iteration formula. Denote this integral by A m ^. 
An integration by parts gives A m ^ = 2m 2 * lk _ 1 An,fc-i- Since A m $ = 2 m+ 2/(2m — 1), we thus obtain 



A 



(2m - l)(2m - 3) • • • (2m - 2k - 1) 
which gives 



(3.18) 



d k U (£,l). 2 k k\m(m - 1) ■ ■ • (m - fc) 



(i£ fc (2m- l)(2m-3)---(2m-2Jfe-l) 
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Therefore 

d k p(0) _ d fe+1 C/ (f,l)| .d fc C/ (£,l) 2 k (k + l)\m(m-l)---(m-k) 

~ d^+ l l?=0 ~ ^ 3 ^ = ° ~ (2m - l)(2m - 3) • • • (2m - 2fc - 3) > ' 

We now use ( I3.17P to prove the existence and uniqueness of the zero of function k = 

0, 1, • • • , (m — 2). First, it follows from ( 13. 17ft that has an odd number of zeros in < £ < 1, 

counting multiplicities. Second, if has more than one zero, it must have at least three zeros. 

Consequently, p^ k+1 '(£) will have more than one zero; so it must also have at least three zeros. 
Repeating this argument, we see that p( m_2 ) must have at least three zeros. This is an impossibility 
since is a polynomial of degree m — 1 because Uq(£, 1) is so. Therefore, has one and 

only one zero for £ e (0, 1) when k — 0, 1, • • ■ , (m — 2). In particular, the = case proves Lemma 
1331 □ 

Lemma 3.4. Equation $3.13\) has a unique solution satisfying U2 = u 3 . The solution is U2 = u 3 = u* . 
The rest of equations $3. II) at the trailing edge are U\ — 1 and x/t = ^1^(1, u*,u*). 

Having located the trailing edge, we now solve equations (13.11) in the neighborhood of the trailing 
edge. We first consider equation ( 13.131) . We use ( 12.351) to differentiate F at the trailing edge u\ = 1, 
u 2 = u 3 = u* 

dF(l,u*,u*) _ dF(l,u*,u*) = 1Q d 2 q(l,u*,u*) Pg^n^u*) ^ Q 

du2 du3 du2du3 duldu^ 

(3-19) = -%|[(1 - - C/o(C, 1)W > , 

where in the second equality we have used (12.221) . (13.14ft and identities = |jr = ^ d ^ 2 Q U3 on 



u 2 = u 3 . The inequality is a consequence of Lemma | 

Inequality (I3.19P shows that equation (13.131) or equivalently (13.121) can be inverted to give u 2 as a 
decreasing function of u 3 

(3.20) u 2 = B(u 3 ) 
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in a neighborhood of u 2 = M3 = u*. 

We will extend the solution ( 13.2011 of equation (13.121) by decreasing w 3 in the region < u 3 < u* < 
Mo < 1 as far as possible. We need to evaluate the derivatives 2 a — - — - and 2 a 3 — - on the 



solution of d3HJ». It follows from (T2T29|) . fICTD and f[2T32]) that 
(3.21) 



(3.22) 



~ /4 m) ] sKp 1 (l,u 2 ,u 3 ) 
<9w 2 £ - (1 - s)if 

%2_ - /flH = sKp 2 {l,u 2 ,u 3 ) 
du 3 K — E 



on the solution of (13.11) . 

We first study the two polynomials pi and p 2 of (12.331) . 

Lemma 3.5. For each < U3 < 1, the polynomial pi(l, U2, W3), as a function of u 2 , has only one 
zero in the region < u 2 < 1, counting multiplicities. Furthermore, pi(l,u 2 ,u 3 ) is positive when u 2 
is on the left of this zero and negative when u 2 is on the right. 

For each < u% < 1, the polynomial p 2 (l,u 2 ,us) , as a function of u 3 , has only one zero in the 
region < u 3 < 1, counting multiplicities. Furthermore, p 1 (l,u 2 ,u 3 ) is positive when u 3 is on the 
left of this zero and negative when u 3 is on the right. 

Proof. We will prove the first part of the lemma; the second part follows from ^2(1,^2,^3) = 
Pi(l,u 3 ,u 2 ). 

The proof of the first part is similar to the proof of Lemma 13.31 We will go through it briefly. 
We first have 

d k p t (1,0, -u 3 ) Q d k p l {l,l,u 3 ) ^ Q 
du\ du\ 

for k — 0, 1,2, ■ • ■ ,m — 2. The second inequality immediately follows from formula (12.331) for p 1 . The 
first inequality is derived from a formula similar to (I3.18p . 

The rest of the proof is the same as the proof of Lemma 13.31 □ 
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We now continue to extend the solution (I3.20p of equation (13 . 12j) in the region < u 3 < u* < u 2 < 1 
as far as possible. When u 2 and u 3 are close to u*, because of (13391) . we have ^ 2 du ^ 3 ' > and 

p., (m) _ 

KH ' > on the solution of (l3~12l) . These along with (l3T2T|) and (E22) show that 
(3.23) p 1 (l,u 2 ,u 3 ) < , p 2 (l,u 2 ,u 3 ) > 

when u 2 and w 3 are close to u*. 

Lemma 3.6. Inequalities Ii3. 23\) hold on the solution of equation li3.1S\) as long as < u 3 < u* < 
u 2 < 1. 

Proof. We first prove the first inequality of (13. 23ft . In view of (I2.29p . the solution of equation (13.121) 
is also governed by an equivalent equation M(l, u 2 , w 3 ) = 0. For each < u% < u*, we study the zero 
of the «2-variable function M(l,u 2 ,us) in the interval u-$ < u 2 < 1. We obtain from formula (12.301) 
for M that M(l, u 3 , m 3 ) = and M(l, 1, m 3 ) > 0. M(l, u 2 , m 3 ) is a decreasing function of u 2 when u 2 
is on the immediate right of u 3 . To see this, we note that pi(l,u 3 ,u 3 ) > for w 3 < u* according to 
Lemma 13.31 This and (12.311) prove that M is decreasing for u 2 on the immediate right of w 3 because 
E — K < for s > 0; so M(1,m 2 ,m 3 ) < for such u 2 . Therefore, M(l,w 2 ,u 3 ) has a w 2 -zero in the 
interval w 3 < u 2 < 1 when < w 3 < u*. Because of the uniqueness of the M 2 -zero of p(l,u 2 ,us) 
according to Lemma |3.5[ we conclude that M(1,m 2 ,m 3 ) has only one zero and that this zero is on 
the right of the zero of p\. Hence, the zero of M is exactly given by u 2 = B(u 3 ) and p\ is negative 
on the solution. This proves the first inequality of (13.231) . 

We now prove the second inequality of (I3.23|) by contradiction. Suppose it first fails at u 2 and u 3 , 
where < m 3 < u* < u 2 < 1; i.e., 



(3.24) p 2 (l, B(u 3 ), u 3 ) > when u 3 < u 3 < u* , p 2 (l, B(u 3 ),u 3 ) = when u 3 = u 3 
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It then follows from (I3.2ip and (I3.22p that solution (13.201) of equation f 1 3 . 1 2 j) has a zero derivative at 
u 3 = u 3 ; i.e., 

B'(u 3 ) = . 

Hence, 

d d _ _ _ d 

—p 2 (l,B(u 3 ),u 3 )\ U3=il3 = —p 2 (l,B(u 3 ),u 3 )B'(u 3 ) + —p 2 (l,B(u 3 ),u 3 ) < , 

where the first term vanishes because of B'(u 3 ) = and the second term is negative according to 
Lemma [3. 51 In view of p 2 (l, B(u 3 ),u 3 ) = 0, this implies that p 2 (l, B(u 3 ),u 3 ) < when u 3 is on the 
immediately right of u 3 . That contradicts (I3.24p . This proves the second inequality of (I3.23p . 

□ 

It follows from fl3T2T|) . and Lemma 3.6 that 

(3.25) at4"'-/4"'] >0| %r'-^'] >0 

ou 2 ou 3 
on the solution of (13.121) . Solution (I3.20p of equation (13 . 1 2[) can then be extended as a decreasing 
function of u 3 as long as < u 3 < u* < u 2 < 1. 

There are two possibilities: (1) u 2 touches 1 before or simultaneously as u 3 reaches and (2) u 3 
touches before u 2 reaches 1. 

It follows from (I2TT^|) and fl2T23|) that 

/4 m) (l, l,u 3 ) > /4 m) (l, 1,« 3 ) for < us < 1 . 

This shows that (1) is impossible. Hence, u 3 will touch before u 2 reaches 1. When this happens, 
equation (13.121) becomes equation (13. 4p . 

Lemma 3.7. Equation ( |ff.^| ) has a simple zero in the region < u 2 < 1, counting multiplicities. 
Denoting the zero by u** , then fi 2 m (1, u 2 , 0) — /i 3 m ^(l, u 2} 0) is positive for u 2 > u** and negative for 
u 2 < u** . Furthermore, < u*** < u** < 1 where u*** is the unique zero ofp\(l, u 2 , 0). 
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Proof. We use ( I2T291 and OUT]) to prove the lemma. In equation (l2~29l) . K - E and E - (1 - s)iiT 
are all positive for < s < 1 in view of (12.111) . 

Denoting the unique zero of w 2 , 0) by u***, it then follows from Lemma [331 that pi(l, «2, 0) > 
when < u 2 < u*** and pi(l,u 2 ,0) < when u*** < u 2 < 1. Since M(l,w 2 ,0) of (I2.30j) vanishes 
at «2 = and is positive at u 2 — 1 in view of (I2.5H2.8I) . we conclude from the derivative (I2.3ip that 
M(l, U2, 0) has a simple zero in < u 2 < 1. This zero is exactly u** and the rest of the theorem can 
be proved easily. □ 

Having solved equation (13. 12ft for u 2 as a decreasing function of u% for < u 3 < u*, we turn to 
equations (13.11) . Because of (13.111) and (I3.25p . the third equation of (13.11) gives u% as a decreasing 
function of x/t for a < x/t < (3, where a = /4 (1, u*, u*) and (3 = /4 (1, u**, 0). Consequently, u 2 
is an increasing function of x/t in the same interval. 

Lemma 3.8. The last two equations of ( Iff. 1\) can be inverted to give u 2 and u% as increasing and 
decreasing functions, respectively, of the self- similarity variable x/t in the interval a < x/t < f3. 

We now turn to equations (13.21) . We first use (12.91) . (I2.10p and (12.121) to calculate the derivative of 
/4 m) of frZ25|) 

d/4 m) ,« , , ,d 2 q ld\ 2 dq 

~cW = 2 [A2 - 2(1 + M2 + M3)] M + 2^a^ 

2{u 2 -l)sK d 2 q 2sK E 2 - (1 - s)K 2 dq 

E - (1 - s)K duj + E - (1 - s)K + ~ (E - (1 - s)K) 2 ^du~ 2 

(-) > -j-^Rl-^g-*,, 

where in the inequality we have used (J5 — (1 — s)K) 2 > i? 2 — (1 — s)fC 2 , which is a consequence of 

(EUJ). 

The polynomial in the parenthesis of (13.261) is connected to pi(l, ^2, 0); indeed, 

f oo7\ n m 2m+l <9 2 g(l,w 2 ,0) dq(l,u 2 ,0) 
(3.27) 0) = — — [(1 - « 2 )— ^ ^— ] • 
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This follows from the identity 

{QO Q\ o <9g(l,u 2 ,0) dq(l,u 2 ,0) , n , 

(3.28) 2m = {l-u 2 ) = hm? 1,m 2 ,0 . 

To see this, taking the derivative of (13.281) and using formula f!2.33j) for pi and equations (12.221) for 
q yield fl3T2~TD . 

To prove (I3.28j) . we use the integral formula (I2.24p for q to calculate both sides of the identity. 
The left equals 

J ((1 + v) + (1 - v)u 2 ) m - x (1 - v 2 )- x l 2 dv . 
The right is 

J ((1 + v) + (1 - v)^)™- 1 (1 - z/ 2 )- 1 / 2 ^ . 

Both sides are equal in view of an easy identity 

m j (l + /i) m ^ 1 (l- / u) 1/2 rf/i= ~J (l + / u) m (l -^)- 1/2 d(i . 

We have therefore proved identity (I3.28p . 

By Lemma [3~51 pi(l,u 2 , 0) is negative for u 2 > u***, where u*** is the unique zero of p\. Since 

o (m) 

u*** < u** according to Lemma [3~71 we conclude from (I3.26P and (I3.27P that ^ > on the solution 
of (I3.2p when u 2 > u**. Hence, the second equation of (13.21) can be solved for u 2 as an increasing 
function of x/t as long as u** < u 2 < 1. When u 2 reaches 1, we have 

z/t = 7 = /4 m) (l,l,0) . 

We have therefore proved the following result. 

Lemma 3.9. The second equation of A3. || ) can be inverted to give u 2 as an increasing function of 
x/t in the interval f3 < x/t < 7. 



m-123/2 



7T 



l + Zi)" 1 "^! -fi) 1/2 d^ 



4?n-l25/2 7r 



1 + //) m (l - /z)- 1 /^// 
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We are ready to conclude the proof of Theorem 3.1. 
The Burgers type solutions (13.51) and (13. 6ft are trivial. 

According to Lemma [3781 the last two equations of ( 13. II) determine U2 and U3 as functions of x/t 
in the region a < x/t < (3. By the first part of Lemma [3.2[ the resulting u\, U2 and M3 satisfy the 
Whitham equations ( 11.61) . Furthermore, the boundary conditions (13. 7p and (13. 8ft are satisfied at the 
trailing edge x — at. 

Similarly, by Lemma [3.91 the second equation of (13. 2p determines U2 as a function of x/t in the 
region (3 < x/t < 7 = 4 m . It then follows from the second part of Lemma 13.21 that u\, 112 and 1*3 
of (13.21) satisfy the Whitham equations (11.61) . They also satisfy the boundary conditions (13.91) and 
(I3.10p at the leading edge x — 7 t. 

We have therefore completed the proof of Theorem 13 . 11 

4. The Minimization Problem 

The zero dispersion limit of the solution of the higher order KdV equation (11. 4p with step-like 
initial function ( 11.71) is also determined by a minimization problem with constraints [SI El [H] 

(4.1) Minimize! — - / / log ^— - ^{ri)^U)drid^ + [ [rix - 4 m r] 2m+1 tU(r])dr]} . 
{^>o, vei 1 } 2n J J ij + fl J 

In this section, we will use the self-similar solution of Section 3 to construct the minimizer for 
m > 2. The m = 2 result has already been obtained in [9]. Even in the m = 2 case, the key 
calculations presented here are different from those in [9]. 

We first define a linear operator 

2tt Jo \V + vJ 

The variational conditions are 

(4.2) Lip = xr]- A m t7] 2m+1 where i/j > , 

(4.3) Lip < xr/ — A m tr] 2m+1 where ip = . 
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The constraint for the minimization problem is 



(4.4) 



The minimizer of (14. ip is given explicitly: 



Theorem 4.1. The minimizer of the variational problem OTIP is as follows: 
(1) For x < at, 



ip{r)) 



-xrj + 4 m (2m + i)tr) 



P m (ri 2 l,u*,u*) 
rj 2 —u* 



(2) For at <x < (3t, 



-x V P { V 2 ,l,u 2 ,u 3 )+4: m (2m+l)triP m (v' 2 ,l,u 2 ,u 3 ) 
yj (l-ri 2 )(u2-ri 2 )(u3-ri 2 ) 



)+4 m (2m+l)fr ? P m (r ? 2 ,l 
V 7 ( 1 -V 2 )(v 2 -'u>2)(v 2 -u>3) 

where Pq and P m are defined in (OOP and u<i and w 3 are determined by equations A3. 
(3) For f3t<x< 4 m t, 



< rj < y/u~3 

Ui < 7] < y/ll^ 
U% <T] < 1 , 



< T] < y/lL^ 

-xP (ri 2 ,l,U2 ,0)+4 m (2m+l)tP m (y 2 ,l,u 2 ,0) 



y/ {l-r) 2 )(y 2 -u 2 ) 



«2 < rj < 1 



where u<i is determined by h3.2\) . 
(4 ) Forx> A m t, 



4>(v) = • 



Proof. We extend the function ip defined on [0, 1] to the entire real line by setting i/)(r)) = for rj > 1 
and taking if> to be odd. In this way, the operator L is connected to the Hilbert transform H on the 
real line [5]: 

14{rj) = / Hil){u)dn where H^irf) = -P.V. / -^dfi . 

J0 ' 71 J-oo V ~ f 1 
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We verify case (4) first. Clearly ip{v) = satisfies the constraint (14. 4p . We now check the 
variational conditions f l4.2ti4.3H . Since ip = 0, 



Lip = < xt] — 4 m tr) 



m +rri 2m+l 



where the inequality follows from x > 4 m t and < 77 < 1. Hence, variational conditions 04.2H4.3p are 
satisfied. 

Next we consider case (1). We write ip(t]) as the real part of gi(i]) for real rj, where 

- /ZT[_a^ + 4 m (2m + ij^MiWll 



2m\ I v ' ' rj^-u* 



9l = y/-l( x - 4 m (2m + l)tr] 2m ) + 



The function g\ is analytic in the upper half complex plane Im(rj) > and gi(r]) ~ 0(1 /rf) for 
large \t]\ in view of the expansion (12.11) for P m . Hence, Hif)(rj) = Im[gi(r])] = x — A m (2m + l)ti] 2m on 
< i] < 1, where H is the Hilbert transform [5J. We then have for < i] < 1 

Lip( v ) = [ V Hrj(n)dii = xt]- 4"V m+1 , 
Jo 

which shows that the variational conditions (I4.2p and (I4.3P are satisfied. 
To prove f]4.4j) . we first claim that 

(4.5) - a + 4 m (2m + ^^Mlh^l > , 

Tj — It* 

for < T] < 1. To see this, we use (12. 131) and (12.251) to calculate a = /4 (1, u*, u*) and (12. 161) and 
( I2.2ip to evaluate P m (l,u*,u*). The left hand side of (I4.5|) equals 

0( 2 1V 2 * w 2 -, * « Uof2 , dq{l,u*,u*) dq{l,u*,u*) 
2{T] — 1)(t7 -u )$(r? ,1,m ,w ) + 2(rj - u )[ — + — + — ] 

= 2(rf - l)[C m U{r,\ 1) - C m U(u*, 1)] + 2{jf - u*)C m U(u*, 1) 

(4.6) = 2[(1 - u*)C m U{u\ 1) - (1 - rf)C m U{rf, 1)] , 

where in the first equality we have used (13.141) and the identity (cf. (3.28) of [2]) 

C m U( V 2 ,l)-C m U(u*,l) 



$(7/ , l,u*,u* 



rf — u* 
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In view of ( I3.16p . we have 
d 



[(l- V 2 )C m U( V 2 ,l)] = C m {(l- V 2 ) 



dU{rfX 



U( V \l)}=p 1 (l, v \ v 2 ) . 



d(r?y K ~ ;,~-»~v l ,~n -t« L v- ■, > dr]2 
This derivative is positive when rj 2 < u* and negative when rj 2 > u* because pi is so according to 



Lemma [3.31 Therefore, (1 — r] 2 )C m U(r] 2 , 1) has a maximum at if 



u ; l 



i.e., (l-^ 2 )C m ^ 2 ,l) < 



(1 - u*)C m U(u*, 1) for < 7] 2 < 1. This together with flUD proves Kty . 

It follows from x < at and inequality ( 14. 5 ft that ^ > 0. Hence, the constraint ( 14.41) is verified. 

We now turn to case (2). By Lemma 3.5, the last two equations of ( 13. II) determine u 2 and M3 as 
functions of the self-similarity variable x/t in the interval a < x/t < (3. 

We write ip = Re ((72) for real 7], where 

v^ll-xriPoir] 2 , 1, U2, u 3 ) + 4 m (2m + l)t^P m (r/ 2 , 1, w 2 , u 3 )] 



#2 



-l(x-4 m (2m + l)tr7 



2m\ 



vV ~ 1) (?7 2 - u 2 ){r] 2 - u 3 ) 
The function g 2 is analytic in Im(ri) > and #2(7?) ~ 0(l/rj 2 ) for large |??| in view of the asymptotics 
(12.21) for Pq and P m . Hence, taking the imaginary part of g 2 yields 



We then have 



x - 4 m (2m + l)t7] 2m < r] < 

X - 4 m (2m + l)tr] 2m - l- xP o(v 2 ^2,0)+i m (2m+l)tP rn ( V 2 ,l,u 2 ,0)]r 1 

x - 4 m (2m + l)t^ 2m 



u 3 <7] < y/u 2 

VU2 < V < 1 . 



(4.7) 



xr? - 4 m t?7 2m+1 < 77 < 

xr? - 4 m t r ^ 2m + 1 < 77 < 1 , 



where we have used 



[-xP + 4 m (2m + l)tP m ]n 
a/(1 - /i 2 )(n 2 - /i 2 )(/i 2 - W3) 

which is a consequence of ( 12. 3p for P and P m . 



dfi = 
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To verify (14. 4p . we derive an integral formula for ip. We use (12.281) and (13.11) to calculate 

riP (ri 2 ,l,u 2 (x,t),u 3 (x,t)) 



^/{l-^){u 2 {x,t)-^)(u 3 {x,t)-^) 

ipx{v) = { 



< rj < y/u 3 (x,t) 
y/u 3 (x,t) < rj < y/u 2 {x,t) 



-r/P (r/ 2 , l,n 2 (x,t),0) 
k ^{l-V^-Mx^W-usixJ)) 



Integrating yields 



(4- 



rx(v 2 )t ■ n P a {r t 2 ,l,u 2 {y,t),uz{y,t)) , 



\/u 2 (x,t) < rj < 1 . 



< 7] < y/u 3 (x,t) 



X(ri 2 )t 





V P ( V 2 ,l,u 2 (y,t),0) 



--dy 



y/u 3 (x,t) <T] < y/u 2 (x,t) 

\/u 2 (x,t) < r) < 1 , 



where xiv 2 ) is defined in Figure 1. The polynomial Po of (12.1ft is linear in rj 2 and has a zero for 
u 3 < rf < u 2 because of (12.31) . P must be positive for rf > u 2 and negative for rf < u 3 . This 
combined with (14. 8ft proves ip > 0; so (I4.4p is verified. 

We now continue to verify the variational conditions ( 14. 2ft and ( 14.31) . Again, we use ( 12.28(1 and 
(13.11) to calculate 

d_ I [-xP (/i 2 , 1, u 2 (x, t),u 3 (x, t)) + 4 m (2m + l)tP m (fi 2 , 1, u 2 {x, t),u 3 (x, t))]f i 

dx V ^-^MM)-/'^ 2 -*.*)) 

_ -/xPqQu 2 , 1, tt 2 (ff, t),u 3 (x, t)) 

~ y/(l-fjP)Mx,t)-tA*)(f-i4(y,t)) ' 

Integrating yields 

[-xP (/i 2 , 1, U2(x, t),u 3 (x, t)) + 4 m (2m + l)tP m (/i 2 , 1, u 2 (a:, t), u 3 (a:, *))]// 



fiP (lJ, 2 ,l,u 2 (y,t),u 3 (y,t)) 



--dy 



'x(^)t y/(l ~ /i 2 ) (My, t) - /U 2 )(/i 2 - ul(y, t)) 
for yju 3 (x, t) < /j, < a/ m 2 (x, t). The single integral in (14. 7p can then be written as a double integral. 
After interchanging integrals and using (12. 3p for P , the double integral is simplified as 

fiP (fi 2 ,l,u 2 (y,t),u 3 (y,t)) d 



X (r?)t JyfiM - fx 2 )(u 2 (y,t) -/i 2 )(/i 2 -ul(y,t)) 
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The polynomial Po is linear in /j 2 and has a zero for u^{y, t) < fi 2 < u 2 (y, t). In view of (12. 3p for P , 
we must have 



/•'/ 



/xP (/x 2 , l,u 2 (y,t),u 3 (y,t)) 
s/Mut) Vi 1 ~ V 2 )My,t) -/i 2 )(/x 2 -ul(y,t)) 



dfi < 



for a/ us(y, t) < r] < a/ M2(y, t). Hence, the integral in (14.71) is positive and this verifies the variational 
conditions (i~2} and (TP) . 

We finally consider case (3). By Lemma 3.6, the second equation of (13.21) determines u 2 as an 
increasing function of x/t in the interval f3 < x/t < 4 m . 

We write ip = Re(gs) for real rj, where 

9S = yzr (x _ 4 - (2m + lWn + ^l-.W,i,..,o) + ^+iW,i,^o)l . 

a/(?7 2 - 1)(?? 2 - M 2 ) 

The function ^3 is analytic in Im(r]) > and #3(77) ~ 0(l/r) 2 ) for large |r/| in view of the asymptotics 
(12.21) for P and P m . Hence, taking the imaginary part of g 3 yields 



We then have 



V(l-»? 2 )(«2-»? 2 ) 

x - 4 m (2m + l)tr] 2m ^T 2 < 77 < 1 . 



^ _ 4 m t 2m+ l _ p -s^+4"(2m+l)*P mrf < 

I4(rj) = I v/(i-m 2 )(« 2 -^) 

[ XT] - A m tT] 2m+1 y/ui < n < 1 , 

where we have used 

/"^ -^o(/i 2 , 1, U2, 0) + 4 m (2m + l)tP m (/i 2 , 1, ua, 0) 

/ : au = , 

A, v/(l-/i 2 )(« 2 -/i 2 ) 

which is a consequence of f !2.3j) for P and P m . 

The variational conditions (I4.2ti4.3l) and the constraint (I4.4p can be verified using the method for 
case (2). 

□ 
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